Demonstration of a universal one-way quantum quadratic phase gate 
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We demonstrate a quadratic phase gate for one-way quantum computation in the continuous- 
variable regime. This canonical gate, together with phase-space displacements and Fourier rotations, 
completes the set of universal gates for realizing any single-mode Gaussian transformation such as 
arbitrary squeezing. As opposed to previous implementations of measurement-based squeezers, the 
current gate is fully controlled by the local oscillator phase of the homodyne detector. Verifying this 
controllability, we give an experimental demonstration of the principles of one-way quantum com- 
putation over continuous variables. Moreover, we can observe sub-shot-noise quadrature variances 
in the output states, confirming that nonclassical states are created through cluster computation. 

PACS numbers: 03.67.Lx, 42.50.Dv, 42.50.Ex 
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Measurement-based one-way quantum computa- 
tion using an offline prepared, multi-party entangled 
cluster state, is a conceptually interesting alternative 
to the standard unitary circuit model of quantum 
computation In the cluster-model, universality 

is achieved through different choices of measurement 
bases, while the cluster state remains fixed. Unitary 
gates are effectively applied at each measurement step, 
corresponding to elementary teleportations d, for 
propagating and manipulating a quantum state through 
the cluster. The cluster model also turned out to 
provide new, potentially more efficient approaches to 
the experimental realization of quantum logical gates, 
especially in the quantum optical setting [HI, [y] . 

A translation of the circuit model for quantum com- 
putation over continuous variables (CV) 0, [1| to univer- 
sal cluster computation with CV was given in Ref. 
The canonical, universal gate set for CV is {J73(A),C}, 
where C = {Z{s),U2{n), F ,Cz} with the momentum 
shift operator Z{s) = exp(2is2:), the phase gates Ui{ki) = 
exTp{iKix^), the Fourier transform operator F, and the 
controUed-Z gate Cz = cxp(2zx (g) x) Through 
concatenation, the full set enables one to simulate any 
Hamiltonian in terms of arbitrary polynomials of the po- 
sition X and the momentum p to any precision Q . 

The same set without the cubic gate t/3, i.e., the set C, 
is still universal for realizing any quadratic Hamiltonian, 
that is, the whole group of Gaussian unitary transfor- 
mations, the analogue to the Clifford group for discrete 
variables (DV). In the case of DV, for example, single- 
qubit Clifford transformations are fully covered by the 
Hadamard gate H and the '7r/4'-phase gate UTr/4 act- 



ing upon the qubit Pauli operators as U^^^zj 0,^/4^ 
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-iXZ = — F; full universality for 
single-qubit transformations would then require, in ad- 
dition, the well-known '7r/8'-phase gate [2[, the ana- 
logue to the cubic phase gate C/3 for CV. Focussing on 



CV, the quadratic gate from the universal set C for 
all Gaussian transformations, maps the Weyl-Heisenberg 
displacement operators Z{s) and A"(s) = exp(— 2isp) into 

Ul{n)Z{s)U2{n) = Z{s) 

uI{k)X{s)U2{h) = e"*'^"'x(s)Z(-Ks), (1) 



in analogy to the qubit '7r/4'-phase gate t/7r/4- The effect 
of the phase gate may be more conveniently expressed 
in terms of the generators of the Weyl-Heisenberg group, 
IJ\{k)xIJ2{h) = X, uI{k)pU2{k) —p + kx. 

In quantum optics, it is well-know that there is an 
exact, finite decomposition of any quadratic unitary into 
single-mode squeezers and beam splitters In this 

quantum optical language, the quadratic phase gate U2, 
together with the Fourier transform F, provides single- 
mode squeezing, and the two-mode gate Cz involves 
beam splitting modulo single-mode squeezing. 

In the cluster-based one-way model, the quadratic gate 
can be fully controlled through the local oscillator phase 
of the homodyne detector Here, we experimentally 
demonstrate this controllability, with a fixed, offline two- 
mode cluster state. We show that a large set of squeezing 
transformations can be achieved by means of this one- 
way phase gate; sequential application of the gate would 
lead to universal single-mode Gaussian transformations 
(where changes of the 1st moments in phase space re- 
quire, in addition, p-displacements Z(s), trivially realiz- 
able through a cluster state for CV [9|). 

The output states of our elementary cluster compu- 
tations exhibit sub-shot-noise quadrature variance; thus, 
nonclassical states are created deterministically through 
cluster computation with the degree of nonclassicality 
fully controlled by the measurement apparatus. There- 
fore, our demonstration differs from previous implemen- 
tations of universal offline squeezing in which dif- 
ferent squeezing transformations require different beam 
splitter transformations to achieve universality. 
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The elementary teleportation step for the case 
of CV [§] is described as follows. First, in the 
ideal scheme (Fig. 1(a)), an arbitrary input state is 
coupled to a single-mode, infinitely squeezed state 
(a position eigenstate |a; = 0)), ?/QND|'0)in|a; = 0)a. 
This results in e^^"'»PA J dx'>p{x)\x)in J dp\p)A/V^ = 
J da;'0(a;)|a;)in|a;)A, where the subscripts 'in' and 'A' de- 
note the input and ancilla modes, respectively. Up to 
local Fourier rotations, the resulting state corresponds 
to a perfect two-mode cluster state, already carrying the 
quantum information to be processed through the cluster 
(i.e., the quantum state |V')in)- 

Next, we measure the observable W {x)pU{x) of mode 
1, where U{x) = exp[i/(i)] is diagonal in the po- 
sition basis and p is the conjugate momentum to x 
{[x,p] — i/2). The quantum state after the measurement 
with outcome po is •\/7rin(Po|C^(^in) / '4'ix)\^)in\x)Adx — 
y/n f in{po\x)inU{x)ilj{x)\x)Adx Z(-po)U{xa)\iP)a- 
After correcting the displacement Z{~pq), we obtain the 
desired state U{x)\ip) in the ancilla mode. Through this 
scheme, in principle, we can apply an arbitrary unitary 
operator U{x) to \ip)in', for nonlinear gates such as the 
cubic gate U3, however, this would require measuring a 
nonlinear observable. Here, we consider detection of the 
whole range of rotated quadratures (all linear combina- 
tions of X and p) , effectively applying the quadratic phase 
gate U2{k) = exp (ikx^) to |V')in, up to a phase-space dis- 
placement depending on the measurement result pq. 

In our optical realization, x and p are quadrature op- 
erators, for the mode operator d = x + ip. The quadratic 
gate U2{n) corresponds to a sequence of rotation, squeez- 
ing, and rotation 11], with Xout = a^in and Pout — 
Pin -l- K^in. Thus, the required measurement corresponds 
to measuring 19,] p + nx — y/l + tan^ 9 {p cos 9 + x sin 9) 
with K — tan^. Using homodyne detection and setting 
the phase of the local oscillator (LO) to 9, we can mea- 
sure (p cos 9 + x sin 9) . Appropriate electric amplification 
of the homodyne results with gain (1 -I- tan^6')^/^ leads 
to the desired measurement of U'^pU . We show this for 
several values of k,: 0, ±1.0, ±1.5, ±2.0, with coherent- 
state inputs. The corresponding LO phases are 0° , ±45°, 
±56.3°, and ±63.4°, respectively. 

In our optical demonstration, we use three squeezed- 
vacuum ancillae. One ancilla is coupled to the input via a 
QND gate (denoted by subscript A) . The QND gate itself 
requires two additional squeezed vacuum states (denoted 
by subscripts B, C). For the QND gate, we employ the 
scheme of Refs. [l3,[l3]. The full input-output relations 
of the scheme including finite-squeezing resources are 
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(a)Schematic of elementary one-way quantum 
computation. 
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(b)Our implementation of one-way quantum computation. 

FIG. 1: Schematic of a one-way quantum gate and our experi- 
mental setup. OPO: optical parametric oscillator, LO: optical 
local oscillator, and EOM; electro-optic modulator. 



Even with the excess noise from the finite squeezing of the 
ancillae, we are able to observe sub-shot-noise quadrature 
squeezing for sufficiently large k. In the remainder of 
the letter, we shall describe the experimental details and 
present the results of the experiment. 

Experimental setup. — A schematic of the experimental 
setup is illustrated in Fig. 1(b). The original source of 
light is a continuous wave (CW) Ti:sapphire laser, whose 
output is 860 nm in wavelength and 1.5 W in power. 
Quantum states at the 1.34 MHz sideband are used in 
our demonstration. 

The experimental setup consists of the following parts: 
preparation of the input and ancilla states, the QND cou- 
pling gate, measurement, feedforward, and, finally, the 
verification measurement. 

The input state, a coherent state at the 1.34 MHz side- 
band, is generated by modulating a weak laser beam of 
about 10 /iW using electro-optic modulators (EOMs). 
We prepare three types of coherent states | 
a = ipin, and a = 0, corresponding to phase modulation, 
amplitude modulation, and no modulation of the laser 
beam, respectively. 

In order to prepare the ancilla states, there are three 
sub-threshold optical parametric oscillators (OPOs), 
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each generating a single-mode squeezed state, whose 
squeezing level is -4.3dB, -4.9dB, and -5.2dB. An OPO 
is a bow-tie shaped cavity of 500 mm in length, contain- 
ing a PPKTP crystal ^16'| . The second harmonic (430 nm 
in wavelength) of Ti:sapphire output is divided into three 
beams in order to pump the OPOs. 

The QND gate basically consists of a Mach-Zehnder 
interferometer with a single-mode squeezing gate in each 
arm [l5| . Each single- mode squeezing gate contains a 
squeezed vacuum ancilla, homodyne detection, and feed- 
forward [H,!!!. The variable beam splitters in the QND 
gate are composed of two polarizing beam splitters and a 
half-wave plate. We can eliminate the QND gate and just 
measure the input states by setting the transmittances of 
the variable beam splitters to unity. At each beam split- 
ter, we lock the relative phase of the two input beams by 
means of active feedback to a piezoelectric transducer. 
For this purpose, two modulation sidebands of 154 kHz 
and 107 kHz are used as phase references. For the homo- 
dyne detection, the LO phase is adjusted in accordance 
to the desired k value; the feedforward displacement is 
carried out with the right gain depending on k. 

To verify the output state, we employ another homo- 
dyne detection. As is well known from optical homodyne 
tomography, we can reconstruct the quantum state from 
the marginal distributions for various phases [l3|- We 
slowly scan through the LO phase and perform a series 
of homodyne measurements. The 1.34 MHz component 
of the homodyne signal is extracted by means of lock- 
in detection: it is mixed with a reference signal and then 
sent through a 30 kHz low pass filter. Finally, it is analog- 
to-digital converted where the sampling rate is 300,000 
samples per second. 

The powers of the LOs are about 3 mW. The detector's 
quantum efficiencies are greater than 99%, the interfer- 
ence visibilities to the LOs are on average 98%, and the 
dark noise of each homodyne detector is about 17 dB 
below the optical shot noise level produced by the LO. 
Propagation losses of our whole setup are about 7%. 

Experimental results. — As mentioned earlier, we carry 
out the experiment with three types of input coherent 
states \a): a — xia {x-m — 1.4), a — ipm {pin — 1-3), and 
a = 0. For each input state, we demonstrate the gate for 
seven different k values: 0, ±1.0, ±1.5, and ±2.0. 

Fig. 2 shows the raw data of marginal distributions 
and the Wigner functions reconstructed via maximum- 
likelihood method [l3| . We show the results for the input 
state with the amplitude in x as an example. Each scan 
contains about 80,000 data points which are uniformly 
distributed in phase from to 27r, and every 20 points are 
plotted in the figure (about 4,000 data points). For k = 
(Fig. 2(b)), the input state is regenerated at the output 
except for some excess noise. For nonzero k (Fig. 2(c, 
d)), we can see that the distribution of the p variable is 
shifted proportional to x, with a proportionality factor k. 
As a result, the output states are squeezed and rotated. 
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(a)Input coherent state. 
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(b)Output for ft = 0. 
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(c)Output for K = 1.0. 




n 271 -4 -2 

phase (j> [rad] 

(d)Output for K = 2.0. 



FIG. 2: Input and output states with several k. Left figures 
show raw data of marginal distributions and right ones show 
the Wigner functions, reconstructed via maximum-likelihood 
method [13 . 



In Fig. 3, the elliptic output Wigner functions for 
K = 0, ±1.0, ±2.0 are shown, where the position, size, 
and shape of each ellipse correspond to the averaged am- 
plitudes and variances. Fig. 3(a, b) are for the case of 
a = a;i„: (a) experimental results and (b) theoretical, 
ideal operation. They agree well in positions and incli- 
nations of ellipses, although the ellipses in Fig. 3(a) are 
thermalized because of the finite squeezing of the an- 
cilla states. We estimate the experimentally obtained 
K via Kact = (pout)/(2;in), and the values obtained are 
Ka^t = 0.00, 0.95, -1.04, 1.94, and -2.02 for theoretical 
values Kth = 0, ±1.0, and ±2.0, respectively. The dif- 
ferences in inclinations between experimental and ideal 
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(a)o = (xin = 1.4). (t>)a = (ideal operation). 
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{d)a = ipin (pin = 1.3). 



FIG. 3: (color). Input coherent state (black circle) and output 
states for several k. We assume Gaussian distributions, and 
show averaged amplitudes, variances, (a, c, d): Experimental 
results for three types of input coherent state |q), where a 
is the complex amplitude (a = a'^' + a), (b): Theoretical 
prediction with infinite squeezed ancillae for the same input 
state as (a). 



Wigner functions are less than 3°. The experimental re- 
sults for the other input states are shown in Fig. 3(c, 
d) . The change of the amplitude in the input states only 
affects the positions of the ellipses; the shapes and incli- 
nations of the ellipses remain the same. We can see in 
Fig 3(d) that the input amplitude in the p quadrature 
(pin) is simply reproduced at the output and is otherwise 
not affected for any k. All of these results are in good 
agreement with the theoretical input-output relations. 

In Fig. 4(a), the fidelities of the experimental out- 
put states compared to the ideal pure output states (i.e. 
without excess noise) are plotted. The fidelity quantifies 
the overlap between two quantum states, and it can be 
calculated as in (V' I Pout IV') in- In the case of infinitely 
squeezed ancillae, unit fidelity is achieved. In the ex- 
periment, excess noises due to finitely squeezed ancillae 
lead to non-unit fidelities. Without quantum resources 
(i.e., using vacuum states for ancillary inputs), the ex- 
perimental fidelity is 0.62 ± 0.01 for k = 0, which agrees 
with the theoretical result 0.63 derived from Eq. 
With squeezed- vacuum ancillae, the experimental fidelity 
is 0.81 ± 0.01 for k = 0, which is much better than the 
case without nonclassical resources. For nonzero k, the 
fidelities decrease as |k| increases, because the squeezing 
of the ideal output state grows compared to that used in 
the ancillary states. Experimental results are very close 




(a) Fidelities of output 
states. 



(b)Variances of squeezed 
quadrature of output states. 



FIG. 4: Fidelities of output states and variances of squeezed 
quadrature, (i): experimental results with squeezed ancillae 
and their theoretical curves derived from Eq. ([2}. (ii): ex- 
perimental results with vacuum-state ancillae and their theo- 
retical curves, (iii): theoretical results with infinite squeezing 
ancillae. (iv): shot-noise limit. The vertical axis in (b) is 
normalized to the shot-noise limit. 



to the theoretical curves which are calculated from the 
experimentally obtained squeezing levels of the ancillae. 

In Fig. 4(b), the quadrature squeezing of our setup is 
plotted. Note that the squeezed quadratures are fragile 
and easily degraded by excess noise. In the case of in- 
finitely squeezed ancillae, squeezing is obtained for any 
nonzero k; for k = 0, on the other hand, the variance 
of the input coherent state is preserved. With finitely 
squeezed ancillae, the excess noises are added to the 
variances of the ideal outputs. Without nonclassical re- 
sources, squeezing below the SNL is, of course, not ob- 
tained for any k. In the case of a squeezing level of the 
ancillae below -2.9 dB relative to the SNL, the output 
state is squeezed for sufficiently large \k\. We can ob- 
serve a noise suppression below the SNL by 0.3 ± 0.1 dB 
for K = ±1.0, 0.8±0.1 dB for k = ±1.5, and 1.0±0.1 dB 
for K = ±2.0. 

In conclusion, we have experimentally demonstrated 
the canonical quadratic phase gate for CV in a small clus- 
ter computation. The gate is fully controlled by the local 
oscillator phase of the homodyne detector. We demon- 
strated controllability for a set of coherent input states 
and we observed sub-shot-noise quadrature variances 
in the output states, verifying that our measurement- 
based gate creates nonclassical states. Concatenating 
this scheme would enable one to realize any single-mode 
Gaussian transformation, efficiently applicable to arbi- 
trary input states including non-Gaussian states. 
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